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I INTRODUCTION 


II 

In  conjunction  with  the  U.  S.  Geological  Survey's  untethered,  un- 
manned underwater  inspection  vehicle  study,  pattern  recognition  tech- 
nology is  being  investigated  as  an  adjunct  to  an  acoustic  control  and 
data  link  to  the  surface.  Pattern  recognition  is  terminology  for 
grouping  observations  so  that  only  the  desired  data  is  retrieved  for 
analysis.  Such  grouping  is  desirable  primarily  to  hasten  the  retrieval 
of  wanted  data  and  to  avoid  storing  and  analysing  unwanted  data.  The 
recognition  of  patterns  of  data  can  be  applied  to  optical  or  acousti- 
cal imaging  or  even  to  physical  measurements. 

The  primary  means  of  inspection  used  by  present  day  submersibles 
is  visual  and  photography.  Optical  imaging  is  therefore  an  obvious 
starting  point  for  recognition  studies.  However,  whereas  the  Navy  is 
studying  optical  imaging  for  the  mid-waters  of  the  ocean,  where  tur- 
bidity or  water  clarity  is  not  a major  factor,  offshore  structures 
are  located  around  the  coastline  of  the  United  States  in  relatively 
shallow  waters  which  have  great  turbidity  variance.  In  order  to  be 
able  to  perform  remote  optical  surveillance  of  underwater  offshore 
oil  drilling  structures  and  associated  oil  pipelines,  it  is  neces- 
sary, as  a first  step,  to  examine  the  expected  quality  of  the  image 
received  by  an  active  (i.e.  artifically  illuminated)  untethered  robot 
viewing  system. 

For  a robot  viewing  system  that  is  acoustically  interrogated  and 
controlled,  the  amount  of  information  that  can  be  transmitted  by 
normal  imaging  techniques  (e.g.  TV,  etc.)  will  be  severely  limited 
due  to  the  limitation  in  the  data  rate  of  an  acoustic  data  link. 
Consequently,  it  becomes  desirable  to  consider  the  incorporation  of 
a degree  of  "intelligence"  in  the  robot  vehicle  in  order  to  reduce 
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the  amount  of  data  that  must  transmitted  by  the  acoustic  link  to 
the  surface. 

For  high  quality  imagery,  where  the  intervening  medium  does  not 
distort  the  image,  existing  computer  assisted  techniques  are  avail- 
able to  extract  key  features  of  the  image  (rejecting  undesired  image 
features)  for  analysis  by  preprogrammed  computer  to  determine  and 
transmit  the  desired  information. 

As  the  offshore  oil  drilling  structures  will  be  operating  at  va- 
rious depths  (few  hundred  to  a few  thousand  feet),  the  variability 
of  the  ocean  imaging  characteristics  becomes  paramount.  In  addition 
to  the  variability  of  the  water  clarity  which  distorts  the  image 
received  by  the  robot  viewing  system,  the  available  illumination, 
the  type  of  object  viewed,  the  type  of  detector  used  all  contribute 
to  the  final  image  received.  As  the  final  image  will  vary  in  clarity, 
the  ability  of  a recognition  computer  to  perform  the  necessary  analy- 
sis before  transmission  to  the  surface  will  also  depend  on  this  image 
clarity.  Consequently,  the  performance  of  a remote  surveillance  sys- 
tem will  depend  on  the  ability  of  the  "intelligence"  built  into  the 
system  to  compensate  for  this  changing  image  quality. 

This  present  work  has  examined  the  various  factors  that  determine 
the  image  quality  in  order  to  subsequently  determine  the  feasibility 
of  incorporating  the  requisite  "intelligence"  into  the  viewing  robot 
vehicle.  In  figure  1,  we  have  depicted  in  flow  diagram  form  the  var- 
iables that  go  into  the  formation  of  the  final  image  as  received  by 
the  robot  vehicle.  The  legend  on  the  left  of  the  figure  indicates 
those  parameters  that  this  current  work  has  addressed. 

The  the  most  important  parameter  (box  5,  figure  1)  is  the  imaging 
characteristics  of  the  sea  water.  As  a result  of  this  work,  we  have 
been  able  to  characterize  it  in  terms  of  three  well  defined  variables. 
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These  are  <X  , the  attenuation  coefficient;  s,  the  scattering  coeffi- 
cient in  the  foward  direction;  N,  the  log-log  slope  of  the  scatter- 
ing function.  This  method  of  characterization  has  been  developed 
during  the  course  of  this  study  and  represents  a departure  from  the 
conventional  characterization  of  the  sea-water  optical  properties 
(see  section  II  below) . The  other  variables  have  been  treated  in  a 
conventional  manner  except  that  in  the  final  image  (box  12,  figure  1) , 
the  result  is  expressed  in  spatial  frequency  terms  rather  than  image 
resolution  terms.  The  motivation  for  this  method  of  characterization 
of  the  final  image  stems  from  the  ability  to  treat  imaging  systems 
as  linear  filters  of  spatial  frequencies  (i.e.  reciprocal  of  resolu- 
tion element) . This  is  exactly  analogous  to  treating  electrical  wave- 
forms as  a mixture  of  temporal  frequencies  except  in  the  optical  case 
it  is  two  demensional.  This  method  of  image  characterization  (Fourier 
Analysis)  simplifies  the  image  analysis  immeasureably  by  obviating 
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the  necessity  of  solving  the  more  difficult  diffraction  integrals 
used  for  describing  the  resultant  image. 

In  the  following  sections,  each  box  of  figure  1 is  addressed  se- 
perately  except  box  13.  This  will  require  the  expected  scenario  and 
optical  configuration  to  be  used  for  viewing  before  the  pattern  re- 
cognition hierarchy  can  be  evoled. 

In  section  II,  we  treat  the  most  difficult  variable,  the  charac- 
terization of  the  optical  properties  of  the  sea-water  in  terms  of 
the  three  simple  parameters  mentioned  above.  In  sections  III  and  IV, 
we  examine  the  degradation  that  may  be  caused  by  the  type  of  detect- 
or used  and  the  mode  of  recording.  A summary  is  presented  in  section 
V of  the  results  of  this  analysis  and  the  the  necessary  steps  to 
be  taken  in  future  work  before  the  feasibility  of  using  a pattern 
recognition  technique  for  underwater  surveillance  can  be  evaluated. 
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II  DERIVATION  OF  SEAWATER  TRANSFER  FUNCTION 


As  indicated  in  the  introduction,  the  characterization  of  the 
sea-water  optical  properties  is  the  most  difficult  variable  to 
quantify.  One  approach  has  been  to  formulate  an  expression  for 
the  spread  function  ( angular  distribution  of  intensity  with  distance) 
for  a collimated  or  focused  beam  in  sea-water.  In  general  no  closed 
closed  form  expressions  has  been  found  except  for  short  ranges  which 
involves  the  integration  of  the  volume  scattering  function.  It  has 
been  found  by  Wells that  the  optical  properties  of  the  sea-water 
can  be  expressed  more  easily  by  working  in  the  Fourier  Domain  (i.e. 
working  with  the  Fourier  transform  of  the  spread  function).  However, 
as  explained  below,  we  have  modified  this  method  of  expression  so 
that  experimentally  derived  scattering  function  can  be  used  directly. 
This  modified  method  allows  the  use  of  the  many  measurements  that 
have  been  made  for  a number  of  years  to  be  used  without  resorting  to 
a mathematical  model  for  the  scattering  function. 

A.  Background: 

Originally,  we  had  intended  to  use  the  closed  expressions  derived 

by  Wells  for  the  sea-water  transfer  function  in  conjunction  with  the 

experimentally  derived  scattering  functions  by  the  Scripps  Institute 
(2) 

of  Oceanography  off  the  coast  of  Southern  California. 

We  have  found  that  the  experimental  data  does  not  correlate  well 
with  the  mathematical  model  for  the  scattering  function  given  by  Wells 
It  became  apparent  that  in  order  to  assess  the  image  quality  in  turbid 
waters,  it  would  be  necessary  to  modify  the  Wells  theory  in  order  to 
utilize  the  experimental  data  directly  without  the  use  of  a mathemati- 
cal model  for  the  scattering  function.  We  have  attempted  to  formulate 
the  modification  using  simple  parameters  resulting  in  the  development 
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of  the  imaging  equations  from  the  first  principles  of  Fourier  Optical 
Theory.  As  a result  of  the  additional  time  that  was  necessary  to  de- 
velop the  imaging  equations  from  ground  zero,  so  to  speak,  we  are  not 
as  far  along  as  originally  planned.  However,  the  equations  that  have 
been  developed  have  general  application  and  only  require  the  experi- 
mentally derived  scattering  function,  for  a given  type  of  ocean,  de- 
tector characteristics  and  object  structure  to  assess  the  image  quality. 

B.  Theory; 

The  underlying  equation  that  governs  imaging  in  incoherent  light 
is  given  by, 


(i) 


where, 

I(x,y)  = resultant  image  intensity  distribution  in  the  image  plane 

CME^)  = object  intensity  distribution  in  the  object  plane 

s(x,y)  = response  of  the  system ( seawater  plus  optics)  to  a 
point  source — the  point  spread  function 

and  where  the  integral  represents  a convolution  of  the  object  distri- 
bution with  the  point  spread  function.  Equation  (1)  can  be  evaluated 
more  readily  by  mathematically  transforming  to  the  spatial  frequency 
plane  by  direct  Fourier  transformation  resulting  in. 


-N_, 


(2) 
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where  Z~  is  the  optical  transfer  function  and  the  represents 

the  Fourier  transform  of  each  distribution.  The  u;'s  represent  the 
conjugate  variables  to  the  coordinates,  termed  the  spatial  frequency 
with  units  of  radians  per  unit  length.  As  can  be  seen,  direct  multipli- 
cation is  less  difficult  to  evaluate  than  the  convolution  integral  of 
equation  (1) . In  addition,  the  spatial  frequency  representation  is  a 
more  meaningful  criterion  to  use  when  evaluationg  resolution  require- 
ments . 

In  underwater  viewing,  it  is  more  convenient  to  cast  the  equations 
in  terms  of  angular  coordinates.  Consequently,  equation  (1)  reads, 


where  the  0's  are  the  coordinates  divided  by  the  range  on  the  object 
side  and  by  the  focal  length  on  th  image  side.  For  example,  the  0* 
would  be  equal  to  x/F.  In  addition,  the  intensities  are  now  expressed 
as  energy  or  power/steradian . Similarly,  equation  (2)  becomes. 


I C coqi 


(4) 


where  theCO's  are  given  interms  of  radians  ( frequency ) /radians  (angle)  . 
1.  Derivation  of  the  Point  Spread  Function 

In  order  to  make  the  problem  of  deriving  the  point  spread  function 
in  seawater  tractable,  we  reverse  the  order  of  equation  (3)  and  con- 
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sider  what  happens  when  a point  source  in  the  image  plane  is  pro- 


jected onto  the  object  plane.  The  geometry  is  depicted  in  Figure  2. 


Consider  an  element  of  area  dA  in  the  x,y  plane,  dxdy.  In  angular 


coordinates,  dA  is  d0  d 0 F . A radiating  surface  in  the  x,y  plane 

x y 


can  be  written  as  noy0(^x/  @y)  where  nQ  is  the  energy/power  density 

radiating  from  the  surface  and  P ( ev,ev  ) is  the  diffuse  reflectivity 

y x y 


of  the  surface  at  the  coordinates  0^,©^.  The  brightness  of  the  sur- 


face becomes  nQyO  ( 0 x, 0y)  assuming  a Lambertian  surface.  Consequently, 


the  amount  of  energy  emanating  from  the  elemental  area  is, 


nc  D(9*,0y)' 


2 2 

The  lens  receives TTD  /4F  of  the  energy  where  D is  the  diameter  of  the 


lens.  Hence  the  energy  received  by  the  lens  is, 


P gq-  riop(6*,Qy\  y F 

L A cT*  TT 


£ _ D^nc.yOOx^y')  clOxdQ) 

L ^ 


Now  consider  a single  ray  along  the  axis  of  the  lens,  represented 


by  ET dvdw/AT (AT , the  area  of  the  lens),  to  be  propagating  toward  the 

L Li  L) 


£ , tj  plane.  Experimentally,  a narrow  collimated  beam  traversing  a 
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length  A in  water  scatters  according  to  O' ( 0 ) , the  volume  scattering 
function.  Here  the  A is  assumed  to  be  small  enough  that  only  single 
scattering  occurs.  At  the  end  of  the  path  length  A.,  the  angular 
distribution  for  the  single  ray  is  0"(  0)  ^ (1  -CxA)  where  the  OC 
is  the  seawater  attenuation  coefficient.  At  the  end  of  the  second  A , 
we  can  consider  the  resultant  distribution  to  be  a convolution  of  the 
angular  distribution  at  the  end  of  one  A with  itself,  i.e.. 


cr(0)  a(i-o<4  0cr(e) a(i  - o<a)  (7) 


where  the  symbol  © means  convolution.  In  this  derivation,  as  experimen- 
tally derived,  the  Cf(0)  has  cylindrical  symmetry  and  hence  a function 
of  0 only.  Continuing  in  this  fashion  until  the  range  R is  reached 
we  have  that. 


5(0}  = (i  — ^C0)AOCT(©)  A. • • . (8) 


where  the  ...  means  that  the  convolution  continues  R/Atimes.  As  the 
energy  of  the  ray  was  E^dvdw/A^,  we  have  that, 


I.,  (el  ^ELs(e)dvaw_ 

al 


* 


- 


where  I is  the  intensity  distribution  of  the  ray  at  the  range  R. 

All  the  other  rays,  since  without  aberration  of  the  seawater  would 
converge  at  the  same  point,  would  produce  the  same  distribution  in 
the  C , XJ  plane  as  the  axial  ray.  Neglecting  the  slight  cosine  de- 
pendence, the  resultant  distribution  is  the  superposition  of  all  the 
rays  and  hence. 


I (0)  = El  s(0) 


Consequently  an  intensity  distribution  in  the  x,y  plane  would  pro- 
duce in  the  £ , ^ plane  the  distribution  given  by. 


Due  to  the  reversibility  of  the  light  rays,  we  can  rewrite  (10)  as. 


1(0*.  S (ii) 


Transforming  (11),  we  have  that, 


X C^Bx  * 00  &y)  ~ ^ f ^ 
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where  the  is  the  transform  for  the  seawater. 

2.  Derivation  of  the 


sw 


From  equation  (8)  and  the  relationship  between  the  C and  the 

sw 

s ( 0 ) we  have  that. 


'Z-CcJq)  ~ F (13) 

where  F { } means  the  Fourier  transform  of.  From  the  convolution 

theorem  we  can  rewrite  (13)  as. 


R/  y ^//\ 

(uj#)  - 0-oi^)  ^(orc ©)4 


(14) 


To  evaluate  (14),  we  use  the  basic  definition  for  the  numerical  con- 
stant e,  i . e.  , 


e = ( 1 + 6) 

6 ->  oo 


>4 


(15) 


With  6 = - for  the  first  factor  in  (14)  and  6 = <S(0)&-  1 for 

the  second  factor,  we  have  that. 


^Cu fc)  - (i  +-  aR?  (1  + €)e+ 1 


.Q"09)R 


€ -»oo 


(16) 
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Consequently,  the  determination  of  comes  down  to  the  determination 
of  the  transform  of  O'  (©). 

3.  Determination  of  the  Transform  of  0(0) 


In  order  to  determine  the  transform  of  (3"  ( © ) , we  usually  must 
assume  some  form  of  the  sacttering  function.  However,  the  experiment- 
al data  indicates  that  the  scattering  function  follows  a log-log 
relationship  with  arbitrary  slope  out  to  at  least  1°.  Consequently, 
we  can  write  that, 

<5'C0')  - K O (18) 


where  K = constant 

N = arbitrary  slope 

Equation  (18)  indicates  a singularity  at  © = 0.  As  experimentally 

o 

the  function  cannot  be  measured  at  angles  much  less  than  . 1 , we 
assume  that  the  log-log  relationship  holds  even  when  © is  less  than 
.1°.  Although  (18)  diverges  at  0=0,  its’  integral  is  finite  as 
will  be  shown  below.  Consequently,  (18)  can  be  used  to  determine  the 
transform  directly. 

Due  to  the  cylindrical  symmetry  of  (3"  (©),  the  two  demensional 
transform  can  be  written  as  a one-  denensional  Fourier-Bessel  transform. 
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g'cs) 


- 2TT 


K / O'NTo(^0de 


where,  J = Bessel  function  of  zero  order 
o 

6.  = limiting  angle  of  integration  (about  1 
Equation  (19)  then  becomes. 


arc©) 


- 2T tK  I w0)d0 


(19) 


(20) 


with  the  restriction  that  N must  be  less  than  2.  If  the  N = 0 or  if  the 
N = 1,  (20)  can  be  evaluated  simply.  We  let  x = u>9  and  substitute 

in  (20)  . Hence, 


Sxe)  = irTK 


CO 


2.-N 


l-fsl 

X X(  x)dx 


(21) 


For  N 


1 


orce) 


N=l 


- 2TTK 

CO 


X,  O ) d X 


(23) 


The  integral  cannot  be  expressed  in  simple  form,  but  because  of  its 

( 4 ) 

importance  it  has  been  tabalated 

The  Schripps  data  indicates  that  N is  not  an  integer  or  a simple 
fraction  but  something  between  1.5  and  1.6.  Consequently,  it  is 
necessary  to  evaluate  (21)  for  the  general  case  of  N<2. 

To  simplify  the  equation  somewhat,  we  let  (1-N)  = n and  rewrite 
(21)  as. 


cs'CeO 


(24) 


We  expand  the  JQ  into  its  series  representation  and  integrate  term 
by  term,  i .e. , 


ste)  = lTr  K 


CO 


n-H 


-Xo 


X"dx 


, Xo 


. Xo 


2X 


X dx.  +-  n"T5  * <*X 
2 ; 


or, 


r>+  * 


~ _ 2.TTK 

&C6)  - — xrr 


n+5 
Xo 


n+t> 

Xo  . 


CO 


Cn  + O 


■FVTTS  "2*.  4* 


(26) 
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Taking  out  x we  obtain. 


'■w  iTTK/^Y+'/ _L_  - J +■  -5 £-.\ 

c T(e)  - ^Tj)  j^n+l)  2 1#Cn  + 3)  2*42  Cn+5)  J 


Substituting  back  for  N,  we  have. 


^-.TTK,  - 


Xo\  J_J ! , . + — — — — - 

cZT/  Ifc-N)  2*  C4-M) 


So  that  the  transform  is  like  J (x,J  with  modified  coefficients 

o o 

depending  on  N,  the  slope  of  the  scattering  curve. 

In  order  to  complete  (28) , we  evaluate  the  constant  K by  considering 
the  integral  of  (18),  i.e.. 


« 3 77"  / (zr'fa)  Qd9  = 2TTK  g'~  d0 


If  we  call  the  integral  of  QT up  to  €■  , S , then. 


. _SU-Nll 
K ITTe4*"0 
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Substituting  back  in  equation  (28)  we  obtain  that, 


(31) 


where  ,2represents  the  series  solution.  Since  = £ , we  finally 

obtain  that, 


<y(e)  = C.2. 


(32) 


This  form  of  the  transform  is  useful  as  the  value  s has  been  calculated 
by  the  Schripps  data.  As  can  be  seen  when  N=  0 and  N = 1,  (32)  becomes 

equations  (22)  and  (23)  respectively. 

Equation  (32)  can  now  be  substituted  back  into  the  equation  for 
equation  ( 17 ) > and  can  be  written  as. 


_ o<  R.  (z-n)S^I  R. 

z-  Cco6)  = e . e 


(33) 


Consequently  £*(  CJp  ) depends  on  CX  , N , €:  and  s,  the  water  characteristi 
and  the  range  R. 

It  is  instructive  to  compare  this  result  with  the  result  given  by 
Wells,  assuming  a mathematical  model  for  Cf  ( © ) . Wells  gives  for  "2^, 


~ . \ _ 5 0 

e . e — 


-e_ 

CU  &o 


— <-o  Qo'j 


where,  0O  is  a characteristic  of  the  scattering  function.  For  OJ  = o , 
2“ = e_olR.esR,  whereas  equation  (33)  also  gives  the  same  result.  As  CO 


goes  to  0,  both  expressions  give 


However,  values  in  between 


do  not  correlate.  For  example,  at  xQ  and  oj  0o  = .5,  the  *2“  of  Wells  gives 
e-c*  R.e*  7^sR  whereas  equation  (33)  becomes  2“  ^ e-**R.e^~  ' SR* 

The  coefficient  of  sR  becomes  .97  for  N = 0 , .98  for  N = 1 and  about 
.99  for  N = 1.5.  As  this  occurs  as  an  exponential,  the  two  expressions 
for  2^  would  deviate  greatly. 

Consequently,  the  Wells  2"  is  a more  rapidly  decaying  function  than 
equation  (33)  and  is  the  result  of  assuming  a value  for  0=0.  Equation 
(33)  represents  a more  gradually  decaying  function  resulting  from  a 
more  peaked  cS"  ( O ) for  small  angles  as  the  experimental  data  indicates. 
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V 


Ill  EVALUATION  OF  DETECTOR  RESPONSE  TO  AERIAL  IMAGE 


In  section  II  we  have  derived  a nearly  closed  form  expression 
for  the  degradation  of  a received  image  due  to  the  sea-water  charac- 
teristics. In  a viewing  system,  this  degraded  aerial  image  must  now 
be  recorded  by  whatever  detector  is  being  used  in  the  system.  How 
the  detector  'reads'  the  aerial  image  will  determine  the  response 
obtained  from  the  overall  system. 

The  expression  that  was  derived  in  equation  (12)  represents  the 
resultant  intensity  in  the  aerial  image  in  spatial  frequency  units. 
We  must  now  consider  what  happens  as  a detector  records  the  image 
to  give  the  system  response.  The  detector  can  be  any  one  of  several 
systems  used  for  detection.  For  example,  the  detector  can  be  a tele- 
vision camera,  a photodiode  imaging  array,  photographic  film,  image 
dissector  tube  or  the  newly  developed  charge  coupled  device.  In 
order  to  make  this  analysis  very  general  so  that  whatever  detector 
is  used,  the  resulting  expressions  will  assess  the  ability  of  each 
type  of  detector  to  reproduce  the  aerial  image,  for  all  detectors 
a recording  medium's  resolution  element  can  be  defined.  For  TV,  it 
is  the  electron  beam  spot  size  and  for  photographic  film  it  is  the 
grain  size  and  so  on.  Due  to  the  finite  size  of  the  recording  ele- 
ment, the  aerial  image  will  be  further  degraded  as  for  each  detector 
there  is  a required  signal  to  noise  ratio  for  the  element  size  ne- 
cessary for  detection. 

Returning  to  equation  (11)  of  section  II  where  the  received  in- 
tensity is  expressed  in  angular  coordinates  of  the  image  plane, we 
wish  to  determine  the  resultant  intensity  expression  in  the  detector 
plane.  We  represent  the  aerial  image  intensity  as  I(  0y,  0y)  and 

the  detector  intensity  as  I ( 0 ,04-).  Let  the  detector  spot  size 

u u 


H 


have  a radius  of ©s which  spreads  out  each  point  of  I(  Qx  ,Qy  ) into 
a diameter  of  2 0$.  This  can  be  expressed  mathematically  as  the 
function  circle ( 0)/TT0s  (i.e.  a uniform  circle  which  intercepts  /jt  Os 

fraction  of  the  incoming  image.  The  angular  coordinates  of  the  re- 
sultant distribution  we  call  I ( e LL  / e*  ) . Hence,  similar  to  the 
object-image  relationship  used  before,  we  have  that, 


TCBu,9t)  = 


VC  Os' 


35) 


Aperture 


As  this  is  a convolution  of  I ( 0x,  0y  ) with  the  recording  spot,  we 
transform  (35)  to  yield. 


I = I 


circle  Cb') 

TT©SX 


(36) 


But  the  circle  ( © ) has  the  transform  TT  ©s  2 0S  . 

Hence  (36)  becomes, 


I C <*>&£.)  - I > ^Oy) 


z T|  Cto  ©s) 
cu  ©5 


(37) 


Equation  (37)  expresses  the  additional  degradation  of  the  image 
spectrum  due  to  the  finite  recording  spot  size.  Due  to  the  manner  in 
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§ 

E 

t 

which  I ( co0K , uJq^  ) was  derived,  equation  (12),  I represents  the  total 
number  of  photons  passed  by  the  system  at  a given  frequency.  For 
example,  at  oj  = o,  I (0,0)  represents  the  total  illumination  or 
total  number  of  photons  in  the  image  plane.  Consequently,  substituting 
the  expression  for  I ( uJgK  , ) , equation  (12)  and  expressing  now 

the  I ( , ^5^)  in  terms  of  photons,  we  have  that, 


N C uJ9u.yCtJ&t}  ~ 


t>5w  pC.^et  .^) 


i&P> 


y 


60  6)  * 


(38) 


In  order  to  calculate  the  resultant  signal  to  noise  photon  ratio 
(S/N) , we  consider  the  effect  of  background  light  in  the  image  plane. 

The  effect  of  background  light  can  be  included  by  treating  the 
background  as  a uniform  object  in  object  space  whose  diffuse  re- 
flectivity is  given  by, 

- TT  rey  Circle  (0)  {39) 


where , 

C5b  = the  scattering  function  value  for  seawater  at  180° 
re^=  the  effective  distance  through  which  backscattering  occurs 
circle(0)  = represents  a uniform  beam  used  for  illuminating  the  object 
It  is  noted  that  since  has  units  of  1/  distance  steradian,  the 
IT  is  necessary  because  of  the  way ^0  was  originally  defined.  The  value 
of  can  be  calculated  by  knowing  the  size  of  the  illuminating 

beam,  the  angle  between  illuminating  beam  and  viewing  system  direction 


-21- 


and  range  R.  Provided  r «R,  it  can  be  shown 


(5) 


that , 


2_  Qmax'K 


(40) 


where  Q^ax.  = angular  radius  of  illuminating  beam  at  the  object 

r = angle  between  illuminating  beam  and  direction  of  viewinq 
The  transform  o f/¥e  ) is  just. 


^8 


- tra’i 


8 


Clirc  L£  C6) 


(41) 


or. 


■O-  _ TT  Ob  Ae  2. 

,42» 

where  Aq  = angular  area  of  the  illuminating  beam  at  the  object.  Since 
the  background  light  traverses  the  same  path  length  as  the  object 
reflected  light,  we  can  write  that. 


^ k°du  >Ufet 


) = 


- n( 


7 USW 


. 2 Xi  Cu>6s) 
toQs 


where,,  reflectivity  spectrum  of  object 


= given  by  equation  (42) 
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(43) 


Equation  (43)  was  derived  in  such  a way  that  N(i.e.  Ng  + Ng)  at 
CO  = 0 represents  the  total  number  of  photons (signal  + background) 
availab) e in  the  entire  image  plane.  This  can  be  seen  by  considering 
the  meaning  of  N when  cl)  = 0.  From  equation  (35),  the  I(  Q 0^)  re- 
presents the  number  of  photons/steradian  as  the  recording  spot  explores 
the  image  plane  with  coordinates  . The  transform  of  this  is  de- 

fined as. 


Apert  one 


Whence)  = 0,  the  transform  is  just  the  integration  of  the  I(  QU,Q^) 
over  the  image  plane. 

In  order  to  define  the  S/N  photon  ratio  properly,  we  recast  equa- 
tion (35)  so  that  I ( Gu. , 0£  ) represents  the  total  number  of  photons 
at  the  coordinates  in  the  recording  plane  rather  than  in  terms 

of  photons/unit  steradian.  This  is  done  by  treating  the  recording  spot 
size  as  just  passing  the  number  of  photons  contained  in  the  record- 
ing spot,  i.e. , 


N (0u,9t)  - //  I y)  circle  (du~O^Ot~^y)dOMd9y  (45) 


Ap«r  -Kn-e 

Equation  (45)  expresses  the  number  of  photons  and  hence  current  or 
voltage  in  the  read-out  system  used  at  the  coordinates  0u.,  0 fc.  Trans- 
forming (45),  we  have  that. 
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w 


N 


(0u,0t)  = TT©s  !(©<,%) 


(46) 


or, 


N (©„©*)  = rr  e?  i(e<,ey)2  t, 

Oo  &s 


(47) 


Equation  (47)  represents  the  resultant  equivalent  number  of  photons 
passed  by  a read-out  system.  It  is  noted  that  it  no  longer  represents 
the  total  number  of  photons  in  the  image  plane. 

It  is  known  that  the  arrival  probability  for  photons  on  a surface 
is  given  by  Poisson  statistics  as. 


■pUiM)  =_yjr. 

xi 


nJ 


(48) 


where,  P(x;N)  = probability  of  receiving  x number  of  photons  when 
N is  the  average  number.  If  the  detector  has  a quantum  efficiency  of 
Q,  it  can  be  shown  that  the  photoelectron  signal  to  noise  ratio  is 


just, 


(6) 


w"V 


(49) 
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Consequently  we  can  write  that, 


("s/aj)  # ('s/a/)  - N Q 


(50) 


In  frequency  space  this  becomes. 


= Q N <51' 


where  the  CD  is  the  auotcorrelation  of  the  transform  of  S/N.  The 
X S/N 

introduction  of  the  autocorrelation  of  the  transform  of  S/N  becomes 
a more  meaningful  criterion  in  terms  of  the  transmission  of  infor- 
mation by  a system  than  speaking  of  a S/N  at  a given  point  in  the 
image.  For  example,  at  CO  = 0,  the  represents, 


(52: 


where  the  brackets  indicate  a time  or  ensemble  average.  Hence  the 

2 

average  of  the  (S/N)  over  the  entire  image  plane  is  just  the  var- 
iance in  S/N.  The  square  root  of  this  function  would  be  the  standard 
deviation  of  S/N  over  the  image  plane.  For  other  spatial  frequencies 
a similar  interpretation  could  be  given. 

From  the  point  of  view  of  eventually  developing  a recognition 
system,  the  ability  of  the  system  to  transmit  frequency  information 
becomes  paramount. 
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In  this  derivation,  we  have  not  attempted  to  include  electronic 
noise  of  the  read-out  system  as  in  general  the  design  of  the  read- 
out system  is  made  so  that  the  shot  noise  is  the  major  source  of 
noise. 
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IV.  EFFECT  OF  DETECTOR  SCANNING 


In  section  III,  we  have  derived  the  expressions  that  relate  the 
effect  of  the  recording  medium's  spot  size  to  the  resultant  response 
from  the  detector.  We  consider  now  the  effect  of  line  scanning  on 
the  final  image. 

Equation  (43)  was  derived  on  the  assumption  that  the  recording 
spot  size  explores  the  image  plane  continously  in  all  directions. 

In  the  case  of  line  scanning  as  in  a television  camera,  (43)  must 
be  modified  to  reflect  this  line  scanning.  We  can  treat  the  scanning 
system  as  if  the  aerial  image  is  sampled,  say  in  the0y  direction.  If 
the  seperation  between  each  scan  line  is  called  0y  (angular  units 
in  the  y direction) , then  it  can  be  shown  using  sampling  theory  that 


equation  (43)  becomes, 

jy/z 

i t 1 

n =^y 


Z J~!  ) evaluated. 

CO  ©S  / N. 


where  is  the  number  of  TV  lines.  Equation  (53)  expresses  the  fact 

that  the  image  spectrum  is  reproduced  at  intervals  of  2TTn/ 0 in  the 

to  Qy  direction  in  frequency  space.  If  the  highest  frequency  of  the 

image  spectrum  to  be  passed  by  the  system  is  COQr  then  it  also  can  be 

shown  that  (53)  goes  over  to  (43)  when  0 = TT/uj  • 

y o 

The  sampling  theorem  can  become  important  when  considering  objects 
with  low  frequency  components  which  are  further  degraded  by  the  sea- 
water transfer  function.  This  will  produce  a lower  frequency  limit 
which  allows  the  use  of  less  scan  lines  per  picture  therby  permitt- 
ing the  available  information  to  be  passed  more  rapidly  than  from 
continous  scanning  (i.e.  continuous  scanning  would  not  produce  add- 
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itional  information  about  the  image) . 

As  an  example  of  how  sampling  theory  permits  information  retrieval 
more  rapidly,  consider  the  object  to  be  a rectangle  (i.e.  simulating 
a cross  beam  that  may  be  present  in  an  underwater  off-shore  oil  dril- 
ling rig) . Additionally  assume  it  is  covered  with  barnacles  giving  a 
diffuse  ref lectivity^  with  an  angular  width  0W  radians  and  an  angu- 
lar length  of  © radians  oriented  along  the  0£  direction,  see  Figure 
1.  The  reflectivity  spectrum  of  the  object  becomes, 

CO  Ql  r\  / 

p.ft,  (S4) 

CO  9c/ 2-  cuS^/2. 


The  highest  frequency  in  the  object  will  be  given  by  the  expression, 


Sin  to 


_ 


= O 


(55) 


or  when. 


6U0  - 


•4TT 


9 


(56) 


If  there  were  no  seawater  degradation,  then  the  image  should  be  samp- 
led at  0„  - w/(,^  = Sw/4  to  retrieve  all  the  information  available 


/ * 


a, 

about  the  object.  However,  because  the  GO  is  multiplied  by  C in 


sw 


equation  (53),  the  highest  frequency  is  now  reduced  according  to. 


3jri  00  0^/2.  - O 


A 
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This  is  now  governed  by  the  ^sw  if  the  limiting  value  of  co  due 

to  ^~sw  is  less  than  the  c<->0*  Consequently,  depending  on  the  values 

of  CK  , s nad  N which  defines  the  Zr  , the  image  need  only  be  sampled 

s w 

at  0 = 'n7 to ( /2r ) • This  value  of  © would  be  greater  than  the  value 

y y 

computed  based  on  the  object  limiting  frequency.  As  no  new  information 
can  be  obtained  by  scanning  closer  than  this,  it  reduces  the  number 
of  scan  lines  necessary  per  picture  or  image. 

Due  to  the  large  number  of  variables  involved  in  assessing  the 
image  quality  and  subsequent  detection  (approximately  15) , it  is 
necessary  to  define  a scenario  of  viewing  and  illumination  configu- 
ration before  the  results  of  this  study  can  be  applied.  Once  the 
image  quality  can  be  assessed  then  the  design  on  the  intelligence 
required  on  a robot  viewing  system  can  be  evolved.  See  figure  1 for 
the  variables  involved  before  recognition  parameters  can  be  designed. 


V.  SUMMARY  AND  RECOMMENDATIONS 


During  this  study,  we  have  established  the  basic  framework  for 
assessing  the  image  quality  and  subsequent  detection  for  underwater 
viewing  systems.  In  particular  we  have  derived: 

a)  the  form  of  the  point  spread  function  in  frequency  space 
(the  transfer  function)  of  sea-water  which  utilizes  available 
experimental  scattering  data. 

b)  the  form  of  the  imaging  equations  in  water. 

c)  the  equations  governing  scattering  background  that  reduces 
the  aerial  iamge  contrast. 

d)  the  equations  governing  the  degradation  of  the  aerial  image 
due  to  the  finite  size  of  the  resolution  element  of  the  re- 
cording system. 

e)  the  frequency  space  equivalent  of  the  S/N  that  results  from 
shot  noise 

f)  the  equations  governing  the  line  scanning  of  the  aerial  image 
It  is  recommended  that  the  next  phase  of  this  study  should  add- 
ress the  following  areas: 

a)  Establish  a viewing  system  scenario 

b)  Assess  the  ability  of  the  following  four  photoelectic  detectors 
to  transmit  image  spatial  frequencies; 

1)  conventional  TV 

2)  photodiode  array 

3)  image  dissector  photomultipier 

4)  charge  coupled  device 

c)  Develop  an  optical  algorithm  tha  can  recognize  preprogrammed 
image  information 


d)  Fabricate  an  optical  simulator  for  sea-water  and  detector 
degradation  of  prescribed  images  in  order  to  assess  the 
optical  algorithms  developed. 

e)  Design  and  fabricate  an  optical  recognition  module  to  be  used 
with  a developed  viewing  system. 
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Reference  In  Text 


1.  Not  considered  in  text 

2.  y)0  (Of,  0^),nuTfcer 

not  considered 

3-^/%(aJee>CL>a^ 

4.  Equation  (42) 

5.  Equation  (33) 

6.  Only  aperture  considered 

7.  D^/4  times 

8.  Equation  (48) 

9.  I (CO  a,  , ) 

10.  Equation  (53) 

11.  Not  considered  in  text 

12.  Equation  (“51) 

13.  Subject  of  future  study 
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12.  RESULTING  RESPONSE 


Image  Frequency  Content 
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Correlation  Spectrum 


13.  PATTERN  RECOGNITION 
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On  This  Spectrum 
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FIGURE  1.  Flow  Diagram  Of  The  Steps  Necessary,  With  The  Variables  Involved,  To  Determine 
The  Image  Quality  With  Subsequent  Detector  Retrieval  In  Frequency  Space. 
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IMAGE  PLANE 


FIGURE  2.  Geometry  Used  To  Derive  The  Imaging  Equations  In  Water  With  Subsequent  Detector  Read-out  (See  Text) . 
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